January 11, 2011 



Measurement of the weak phase a from — > ai(1260)^7r^ 

decays 



SiMONE StRACKtJi] 

Dipartimento di Fisica, Universita degli Studi di Milano 
and INFN Sezione dz Milano - 1-20133 Milano, ITALY 



We present the measurement, performed by the BABAR Collaboration, 
of the weak phase a from the time dependent CP asymmetries in 5° — )■ 
ai(1260)^7r^ decays. The model error induced by penguin contributions 
to the —7- ai(1260)^7r^ channel is estimated from an SU(3) analysis 
of the branching fractions of S — )■ ai(1260)K, B — )• ^1(1270)71, and 
B ii'i(1400)7r decays. 
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1 Introduction 



The measurement of the CKM angle a at present-day 5-factories rehes on the analysis 
of time- dependent CP violating asymmetries in tree-dominated h — )■ mid transitions, 
such as — )• TT+TT^, p~ ^ P^T^^ 1 ai(1260)^7r^ (charge-conjugated reactions are 
implied throughout the text). The extraction of a is limited by the penguin contri- 
butions to the decay amplitude, which shift the value of the phase measured from 
the time distribution of decays by an amount Aa that has to be determined from 
the experiment. One of the strengths of the 5-f actor ies lies in their ability to use 
multiple approaches to the measurement of a, allowing for a better control on model- 
dependent estimates of the penguin contributions by comparison with data in many 
channels. Independent measurements of this angle in different channels also help to 
resolve discrete ambiguities that emerge in the extraction of a. 

The angle a can be measured from CP violating asymmetries in decays of neutral 
B mesons to non-CP eigenstates [T], such as p^v:^ and ai(1260)^7r''^. For the p^vr^ 
final state, a can be extracted without discrete ambiguities by looking at the time- 
dependent asymmetries in different regions of the tt+tt^tt^ Dalitz plot. At the present 
level of statistics, this approach cannot be applied to the four-particle final state 
resulting from B^ — )■ ai(1260)^7r^ decays. Nevertheless, the analysis of the time- 
dependent CP asymmetries allows to derive an effective value aeff, which can be 
related to a under the SU(3) approximate symmetry by measuring the branching 
fractions of a set of auxiliary B decay channels: B — )■ ai(1260)i^', B — > /Ci(1270)7r, 
and B — )■ i^'i(1400)7r [2]. In the following, we report the determination of a in the 
5° -)■ ai(1260)^7r=F decays, with the data collected by the BABAR detector at SLAG. 

2 Branching fraction of — > ai(1260)^7r^ decays 

The —7- ai(1260)^7r^ channel was observed by BABAR in 2006 by reconstruct- 
ing the decay of the ai(1260) axial vector meson (henceforth denoted as ai) into 
the dominant pvr channel. The signal contribution is separated from background by 
means of an unbinned maximum-likelihood (ML) fit to a set of five discriminating 
variables. Two kinematic variables, the energy substituted mass niEs = — p\ 

and the energy difference AE = Eb — a/s/2, where the B four-momentum {Eb,Pb) 
is defined in the e~^e~ center-of-mass (CM) frame, allow to discriminate between cor- 
rectly reconstructed B candidates and fake candidates resulting from random com- 
bination of particles. Topological variables, combined into a Fisher discriminant J-", 
provide further distinction between the jet-like shape of continuum e+e~ — )■ qq events 
(g = u,d,s,c), which is the most abundant source of background, and the more 
isotropic B decays. The two remaining variables characterize the resonant behavior 
of the reconstructed three-particle system in the final state: the 7r^7r^7r+ invariant 
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mass niai and the cosine cos 6h of the angle between the momentum of the bachelor 
pion and the normal to the plane described by the resonant three-pion system, in the 
Oi rest frame. The cos 9h distribution allows to discriminate between different 
hypotheses for the three-pion resonance. The lineshape parameters for the ai meson 
are left free in the fit, to minimize systematic uncertainties. 

A signal yield of 421 ib48(stat.) events is extracted from the fit to the BABAR data 
(218 X 10^ BB pairs), which corresponds to a branching fraction B{B^ — t- afn^) = 
(33.2 ± 3.8(stat.) ± 3.0(syst.)) x 10~^, assuming a 50% branching fraction for the 
ai(1260)^ — )■ 7r"'"7r+7r~ decay [3J. These results are in good agreement with the branch- 
ing fraction extracted by Belle, (29.8 ± 3.2(stat.) ±4.6(syst.)) x 10"^ [4j. 

3 Time-dependence of — > afn^ decays 

With a sample of 384 x 10^ BB pairs, BABAR performed a ML fit to 29300 selected 
events, resulting in a signal yield of 608 ± 53(stat.) (461 ± 46(stat.) events with their 
flavor identified), and measured the time distribution of the B'^ — )■ afrc^ decays 

/;^(At) oc (l±Acp)|l + g[(5± A5)sin(AmrfAt) + (C ± AC) cos(AmdAt)] |, (1) 

where Am^; = 0.502 ± 0.007ps~^ is the B^ — B^ mixing frequency, and q = +1 (— 1) 
if the other B in the event decays as a B^ (B^)- The observed time- dependent rates 
and asymmetry are shown in Fig. [1] (a-c), and take into account the At resolution 
function and the dilution from incorrect flavor assignment. They correspond to Acp = 
-0.07 ± 0.07 ± 0.02, S = 0.37 ± 0.21 ± 0.07, AS = -0.14 ± 0.21 ± 0.06, C = 
-0.10 ±0.15 ±0.09, and AC = 0.26 ±0.15 ±0.07 0, where the first error is statistical 
and the second systematic (dominated by the modeling of the signal distributions and 
by CP violation in the BB background). Linear correlations are at the 0(%) level. 
These parameters can be related to the effective value tteff = « ± Aa, by the relation 

S±AS = ± (C ± AC)2 X sin (2aefr ± 5) , (2) 

where 5 is the strong phase between the tree amplitudes of B^ decays to a^^vr" and 
aiTi^ . The strong phase can be averaged out to yield a^s with an eightfold ambiguity 
in the range [0, 180]°, which can be reduced by assuming 5 ^ 1, as suggested by 
factorization [2]. The selected solutions are a^s = (11 ± 7)° and ctefr = (79 ± 7)°, 
where the error is statistical and systematic combined. 

4 SU(3) analysis and B aiK, B Kitv 

The effect of penguin pollution Aa = — a can be evaluated from auxiliary mea- 
surements by introducing flavor-symmetry arguments. Since an isospin analysis is 
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Figure 1: (a-c): Projections onto At of data (points) for (a) 5°, (b) 5^ tags, and (c) 
the asymmetry between and 5^ tags [5]. (d-f): Continuum-background subtracted 
projections of the data (points) on mKnn for (d,e) S° and (f ) events: (d,f ) events 
with 0.846 < rriKn < 0.946 GeV and (e) events not included in (d,f) with 0.500 < 
mTTTT < 0.800 GeV. The sohd hne is the sum of the fit functions for the decay modes 
/s:i(1270)7r + /s:i(1400)7r (dashed), ir*(1410)7r (dash-dotted), K*{892)'K'k (dotted) ^. 
The dashed curve is normahzed to (d) 545, (e) 245, and (f) 141 events. (g,h): 68% 
(hght) and 90% CL (dark) — (f) regions in (g) 5° and (h) 5+ decays to Kitt [9]. 



not feasible [6], our approach is based on a set of SU(3) relations |2], that allow to 
estimate the size of penguin amplitudes from the branching fractions of the AS* = 1 
partners of the 5° — )■ afn^ decays: B — > aiK and B — )■ Kiati, where the Kia 
state belongs to the same SU(3) octet as the ai meson. This approach is effective 
because in these channels the penguin amplitudes are enhanced by a CKM factor A 
(A ~ 0.23) with respect to the AS" = transitions B — )■ oivr. Bounds on Aa can be 
derived from the following ratios of CP-averaged rates: 

where the ratios of the decay constants parameterize factorizable SU(3) corrections. 
Nonfactorizable contributions to AS* = and AS* = 1 transitions from exchange and 
weak annihilation diagrams, respectively, are neglected ||2j. In the above expressions. 
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/aj and fxiA obtained from the study of r decays |7] . 

With an analysis similar to the one for the — )■ afn^ channel, BABAR measured, 
from a sample of 383 x 10^ BB pairs, the branching fractions of B^ — )• aiK~^ and 
B+ -> atK° decays: B{B° aiK+) = (16.4 ± 3.0(stat.) ± 2.4(syst.)) x 10"^ and 
B{B+ atK°) = (34.8 ± 5.0(stat.) ± 4.4(syst.)) x 10'*^ [Sj. 

The KiA is a mixture of the -ft'i(1270) and -ft'i(1400) axial vector mesons, with a 
mixing angle 9 = 72°: \Kia) = 17^1(1400)) cos6' - |iri(1270)) sin 6*. Both resonances 
decay to Ktttc through similar intermediate resonances and are characterized by over- 
lapping mass distributions, and sizeable interference effects are thus expected. The 
contribution of the Kia state can be isolated by extracting from data the combined 
branching fraction of B decays to ii'i(1400)7r and ii'i(1270)7r, and the relative magni- 
tude (r = tan '(9) and phase (0) of -B — )■ i^i(1270)7r and B — > 7^1(1400)7? amplitudes. 

In the analysis recently performed by BABAR with the final data sample of 454 x 10^ 
BB pairs [9], the combined i^'i(1270) and i^'i(1400) signal is parameterized in terms 
of a two-resonance, six-channel ii'-matrix model [lU] in the P- vector approach [11]: 
the i^'-matrix describes the propagation and decay of the Ki resonances, while the 
P-vector effectively parameterizes the production of the Ki system, along with a 
recoiling bachelor pion, in B decays. The decay couplings and the mass poles are 
determined from the results of the partial wave analysis, performed by the ACC- 
MOR Collaboration, of the diffractively produced Kini system [10]. The production 
parameters are extracted from BABAR data by means of a ML fit to AE, tties, ^ i 
cos 6h, and the invariant mass of the resonant Kittt system {rnKTm), which provides 
sensitivity to the individual contributions of the Ki resonances. 

The continuum-background subtracted rriK-K-K distribution in data is shown in 
Fig. [1] (d-f). Including systematic uncertainties, dominated by the effect of in- 
terference between the Ki and the non-resonant components, the combined signal 
branching fractions are i3(P° Ki{mQ)+Ti- + is:i(1400)+7r-) = 31^? x 10"^ and 
B{B+ is:i(1270)°7r+ + is:i(1400)°7r+) = 29l?? x 10"^. The information about 
the fraction and phase of the two resonances (Fig. [1] (g,h)) is used to calculate 
i3(P° ^ K+aT^-) = Utlo X 10-6 and B{B+ K^^n+) < 36 x 10-^ where the 
latter upper limit is evaluated at the 90% confidence level (CL) [S]. 

5 Bounds on Aa 

The bounds on |Aa| are derived by inverting the relations [2] 



A Monte Carlo method is used to derive the 68% and 90% CL upper limits for 
the bounds: replicas of the input quantities are generated from the experimental 




(5) 
(6) 
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distributions, and for each simulated set of values the above system of inequalities is 
solved. This study yields the bound |Aa| < 11° (13°) at the 68% (90%) CL, and the 
final result a = (79 ± 7 ± 11)° for the solution compatible with the CKM global fits, 
where the first error is statistical and systematic combined and the second is due to 
penguin pollution. 

The presence of non-factorizable SU(3) breaking effects can be tested, e.g., at 
LHCb or at a Super i?-factory, by studying auxiliary decay channels such as K^K^. 
The impact of such corrections can be estimated by writing p'_^ = —c±Xp±, where p'j_ 
(p±) is the penguin amplitude in AS = 1 {AS = 0) transitions, and the departure 
of c-t from 1 quantifies the amount of non-factorizable SU(3) breaking. For c± = 0.7, 
the bounds are expected to increase by about 3°. Finally, a full SU(3) fit may provide 
an experimental test of 5 ^ 1. 

6 Conclusions 

BABAR has measured the CKM angle a from the B'^ — )■ afn^ channel, and has 
obtained a value a = (79±7±11)°. This independent determination of a is consistent 
with the world average of the pp, pvr, and tttt channels and with CKM global fits. 
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